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Abstract 

We review the classical and quantum theory of the Pais-Uhlenbeck oscillator 
Q_i as the toy-model for quantizing f(R) gravity theories. 

^: 

1 Introduction 

It is commonplace in Quantum Field Theory (QFT) that a QFT with higher 
(time) derivatives is believed to be doomed from the point of view of physics, 
because of ghosts or states of negative norm, and thus it should be dismissed. 
The standard reference is the very old result (known in the literature as the 
Ostrogradski theorem pQ) claiming a linear instability in any Hamiltonian 
system associated with the Lagrangian having the higher (ie. more than 
one) time derivative that cannot be eliminated by partial integration. 

The key point of the Ostrogradski method [lj is a canonical quantization 
of the clasically equivalent theory without higher derivatives via considering 
the higher derivatives of the initial coordinates as the independent variables. 

The interest in the higher-derivative QFT was recently revived due to 
some novel developments in the gravitational theory, related to the so-called 
/(i?)- gravity theories - see eg., ref. [2] for a review. The f(R) gravity 
theories are defined by replacing the scalar curvature R in the Einstein 
action by a function f(R). The f(R) gravity theories give the self-consistent 
non-trivial alternative to the standard A-CDM Model of Cosmology, by 
providing the geometrical phenomenological description of inflation in the 
early universe and Dark Energy in the present universe. Despite of the 
apparent presence of the higher derivatives, a classical f(R) gravity theory 
can be free of ghosts and tachyons. A supersymmetric extension of f(R) 
gravity was recently constructed in superspace [3]. 
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Already the simplest model of (R + R 2 ) gravity [5] is known as the viable 
model of chaotic inflation, because it is consistent with the recent WMAP 
measurements of the Cosmic Macrowave Background (CMB) radiation [5]. 
Its supersymmetric extension was recently constructed in refs. [HE]. 

On the one side, any quadratically generated (with respect to the cur- 
vature) quantum theory of gravity has ghosts in its perturbative quantum 
propagator [8J. However, on the other side, any f(R) gravity theory is 
known to be classically equivalent to the scalar-tensor gravity (ie. to the 
usual quintessence) [UJ QUI E], while the stability conditions in the f(R) 
gravity ensure the ghost-and-tachyon-freedoom of the classically equivalent 
quintessence theory [12} I13j. It now appears that in some cases the pres- 
ence of the higher derivatives may be harmless [13] • It also gives rise to the 
non-trivial natural question of how to make sense out of the quantized f(R) 
gravity? 

The f(R) gravity theories are just the particular case of the higher- 
derivative quantum gravity theories which have been investigated in the 
past. They were found to be renormalizable [15] and asymptotically free 
|16| . A generic higher-derivative gravity suffers, however, from the presence 
of ghosts and states of negative norm which apparently spoil those QFT 
from physical applications. However, the issue of ghosts and their physi- 
cal interpretation deserves a more detailed study. The complexity of the 
higher-derivative gravity is the formidable technical obstacle for that. It is, 
therefore, of interest to consider simpler QFT as the toy-models. 

Similar features (like renormalizability and asymptotic freedom) exhibit 
the quantum Non-Linear Sigma-Models with higher derivatives, which have 
striking similarities to the higher-derivative quantum gravity |17[ [T8| I19j . 
However, even those QFT are too complicated because of their high degree 
of non-linearity. 

Perhaps, the simplest toy-model is given by the Pais-Uhlenbeck (PU) 
quantum oscillator in Quantum Mechanics [20]. As was demonstrated by 
Hawking and Hertog [2T], it may be possible to give physical meaning to the 
Euclidean path integral of the PU oscillator, as the set of consistent rules for 
calculation of observables, even when "living with ghosts". The basic idea 
of ref. [21] is to abandom unitarity, while never producing and observing 
negative norm states. 

The idea of Hawking and Hertog found further support in refs. [22} [23] 
where the physical propagator of the PU oscillator was calculated by using 
the van Vleck-Pauli approach (the saddle point method for the Euclidean 
path integral) and Forman's theorem [24] . In this paper we systematically 
review the classical and quantum theory of the PU oscillator from the first 
principles, along the lines of refs. [21 EH 1221 123] . 

Some mathematical connections between the higher-derivative particle 
models, higher spins and noncommutative geometry can be found in ref. [25]. 
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2 Ostrogradski method with higher derivatives 

Consider a one-dimensional mechanical system with the action 

S[q] = J dtL(q,Dq,--- ,D n q) (2.1) 



in terms of the Lagrange function L of q(t) and its time derivatives, where 
A 

dt ■ 



n > 2 and D = 4.. The Euler-Lagrange equation reads 



- r)T 

g(- D %^ = ° (2 - 2) 

The Ostrogradski method [lj gives the Hamiltonian formulation of the 
higher derivative Lagrange formulation by introducing more independent 
variables. 

The independent generalized coordinates Qi are defined by 

Qi = D l - l q (i = l,...,n) (2.3) 

The generalized momentum P n is defined by 

9L ;1 =Pn (2.4) 

D n q=A 



d(D n c 



There are n + 1 independent variables {Qi, ■ ■ ■ , Q n , Pn} that are in corre- 
spondence to the n + 1 variables {D°q, • • • , D n q} of the higher derivative 
action (|27TD . 

By solving eq. (|2.4p with respect to A = D n q (assuming that it is possi- 
ble), one gets 

D n q = A(Q 1: ■■■ ,Q n , P n ) (2.5) 

Therefore, the Lagrange dynamics can be represented in terms of the n + 1 
independent variables {Qi, • • • , Q n , P n } as 

L = L(Qi, ■ ■ ■ ,Q n ,A(Qx,--- ,Q n ,P n )) (2.6) 

A Legendre transformation is used to pass from the Lagrange formulation 
to the Hamiltonian one. With the generalized coordinates {Qi,-- - ,Q n } 
and the generalized momentum P n as the independent variables, the total 
differential of the Lagrangian is given by 

dL 

dQj + P n dA 



dL = E 



Di~lq = Q i 
D n q = A 



dL n dL 

!Tq dQ ' + ^W^) dQ l +P " dA 

j— ^ 



dL ,„ 2-1 dL 



^^W* + g WQi dQ]+l + PndA 

(2.7) 
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where we have used eqs. (|2.2[) and (|2.4p . and 



j=i 



Let us now define the n — 1 generalized momenta as 



^ = E(-^7^ (< = l.-,n-l) (2.9) 



J =4 

They satisfy the relations 

dL 



d(l>q] 

Therefore, eq. (j2.7|) can be rewritten to the form 

r-n—l n 



P t + DP i+1 (2.10) 



d 



Y,Pi{DQ l )+P n A-L =-Y,(DPi)dQi + Y.( D Qi) dP i {211 ) 



~i=l i=l i=l 

Equation (|2.1ip gives rise to the Hamiltonian in the form 

n-l 



H = Y J Pj{DQ j ) + P n A-L (2.12) 



The Hamilton equations of motion are given by 

dH dH 

DQ, = — and DP, = -— (2.13) 

3 PU oscillator 

The PU oscillator [2D] is an extension of the harmonic oscillator with the 
higher time derivatives, and is the particular case of the higher-derivative 
theory introduced in Sec. 2. The special features of the PU opscilator are 

(i) the equation of motion is linear. 

F(D)q = (3.1) 
where F is a linear differential operator; 

(ii) the F is polynomial (with respect to D) with constant coefficients: 

n 

F{D) = Y,caD l (3.2) 

8=0 
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where ao, ■ ■ ■ , a n are the real constants; 

(hi) there is the time reversal invariance with respect to t — > —t. Hence, the 
polynomial F has only even powers of the time derivative D. 
The Lagrangian of the one-dimensional PU oscillator reads 



L(q, Dq,-.., D n q) = - £ |(Z^) 2 (a + 0, a n ± 0) (3.3) 



i=0 



where aj (i = 0, ■ ■ ■ , n) are real constants. The Euler-Lagrange equation of 
motion is given by 



=£(-£>)« 



i=0 
- ao 



-aiD l q 



^0 a ° ■ 



(3.4) 



Accordingly, the differential operator F(D) reads 



2i 



F{D) = Y J {-^D 

i=0 



The equation of motion can be rewritten to the form 

F(D)q = 

The PU Lagrangian takes the form (up to a boundary term) 

L = -^qF(D)q 



(3.5) 



(3.6) 



(3.7) 



The differential operator F(D) can be brought to the factorized form 

(3.8) 



where the constants ooi (i = 1,- • • , n) are the solutions (roots) of the equa- 
tion F(iuj) = 0. Let us introduce n new operators 



n ( n 2 \ 
G^) = n i + ^r 

j=i V 3 J 



(i = ,n) 



(3.9) 



and define the n generalized coordinates as 

Q i = G i (D)q (i = l,---,n) 



(3.10) 
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Those generalized coordinates Qj are called harmonic coordinates. By using 
the harmonic coordinates, the PU Euler-Lagrange eq. (|3.6|) can be rewitten 
to the n equations 

D 2 ' 



1 + 



Qi = 



(3.11) 



It means that the PU oscillator can be interpreted as n harmonic oscillators. 
Accordingly, the PU Lagrangian (|3.7p can be rewritten to the form 



where the n constants rji have been introduced as 



Vi = w 



dF 



d(D 2 ) 



D 2 -- 



To prove eq. (|3.13p . we first notice that it amounts to 

n 

J>Gi(£>) = l 

i=l 

By the definiton of G{D) in eq. (|3.9p we have 



G l (D 2 = -co]) 



n 

k^i 

n 



0,7 



k=l 



LCI 



so that 



i2 



i= i k=i \ 



indeed. Therefore, the constants ^ are given by 



n>-3 



Next, we prove that 



dD 2 



D 2 = 



fc=i 

kjii 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



(3.16) 
(3.17) 



(3.18) 



(3.19) 
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By the use of eq. (|3.8p we find 



CLF _ d yr ( 

dD 2 ~ dD 2 + u y 
j=i j 



A 1 A, d\ 

fc=i j=i j 



n 

Y^-G^D) (3.20) 



so that 



dF 



dD< 



D' 2 -- 



3=1 J 

7=1 j fc=l ~fc 

k^i 

i n 2 

(,1. X X x 1,1. 



k^i 

Equation (|3.19p is now confirmed and, hence, via eq. (|3.18p also eq. (|3,13l) 
follows. 

In terms of the harmonic coordinates (|3.10j) . the Lagrangian L, 

2 U 

with the constants r/j given by eq. (|3.13j) . can be rewritten to the form 

n 

£ = TE"«bW 2 -^) (3-23) 

1=1 1 

up to a boundary term. 

The Lagrangian (|3.23p is just a sum of the Lagrangians of n harmonic 
oscillators. Hence, similarly to a free system of n particles, we can change 
the Lagrangian formulation into the Hamiltonian formulation. We define 
the generalized momenta Pi by taking the harmonic coordinates Qi and the 
velocities DQi as the Lagrange variables, 

dL 



d(DQi) 

2 DQi (i = l,---,n) (3.24) 



aoVi 
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The system of n free particles does not have higher derivatives, so its 
Hamiltonian is 

n 

H = Y J PiiPQ i ) - L (3.25) 

i=l 

Equations (1333]) and (ET2l|) imply 

By rescaling the harmonic coordinates and the generalized momenta as 



Ql ^Q l = ^M Qi and Pl -> P t = ^Mp (3.27) 
Ui r]iV a o 

we get the final Hamiltonian 

2 feM 

The presence of both positive and negative values of the constants r/j in 
the Hamiltonian implies both positive and negative values of energy. The 
constants rji are given by eq. (|3.18j) . If satisfy i < j => < £jj, the 
constants r/j are positive for the odd number i, and are negative for the even 
number i. Therefore, the Hamiltonian is 

H =\jZ{-^- 1 (Pr 2 +^Qr 2 ) (3-29) 

This Hamiltonian can be interpreted as that of n harmonic oscillators, with 
the positive and negative energy levels appearing alternatively. Because of 
that reason, the PU oscillator has an instability (for any interaction). It is 
related to a possible ghost state of negative norm in PU quantum theory 
(see Sec. 6). In what follows we consider the simplest case of PU oscillator 
with n = 2 only. 



4 PU oscillator for n = 2: explicit results 

Let us consider the Lagrangian 

H(!) 2 -^>-t(§) 2 h— *») <«> 

with a scalar potential V(q). In the case of the PU oscillator, the potential 
V(q) is a quadratic function of q. Since the (mass) dimension of time is — 1 



8 



(in the natural units K = c = 1), the dimension of the Lagrangian L is 1, 
the dimension of q is —1/2, and that of the constant a is —1. 

Let the trajectory q be a sum of the classical trajectory q c [ and the 
displacement q, ie. q = q c i + q, where the classical trajectory q c i is a solution 
to the equation of motion (EOM) with the boundary conditions [21 J 

*/: q(0)=q , q(T) = q T , q(0) = qo, q(T) = q T (4.2) 

where the dots above stand for the time derivatives. 

With the boundary conditions (|4.2p . the boundary condition of q is 

J: g(0) = 0, q(T) = 0, £(0) = 0, ?(T) = (4.3) 

The action of q c [ + g is given by 

S[fe+g] = S[q cl }+ J* dt^-V(q cl +q) + V(q cl )+qT\q cl )-YQ 2 ) (4-4) 
where we have introduced the notation 

(4.5) 

q=q c i 

In eq. (|4.4j) the term V(q c i + q) — V(q c i) — qV'{q c i) represents the gap 
between the full action S[q] and the classical action S[qd], which generically 
depends on both the classical trajectory q c i and the displacement q. After 
expanding the scalar potential V in Taylor series, 

V(q d + q) = V(q cl ) + qV(q cl ) + ^q 2 V"{q cl ) + ■■■ (4.6) 

we find that, when the second derivative V" is constant, the gap V(q c i + q) — 
V(qd) — qV'(qd) does not depend on the classical trajectory q c \. It is the 
case when the potential V is a quadratic function of q, like the PU oscillator. 

In the path integral quantization (sec. 7), the gap between the full ac- 
tion and the classical action is a quantum effect. When the potential is 
a quadratic function (like that of the PU oscillator), that quantum effect 
does depend on q, but does not depend on the classical trajectory. In what 
follows, we only consider a quadratic function for the scalar potential in the 
form 

V(q) = ^q 2 (4.7) 

ie. the scalar potential of a harmonic oscillator with the mass m > 0, The 
Lagrangian is given by 

Lpu = \q 2 - ^-q 2 - (4.8) 



rr/f s dV 

V {Qd) = -dq- 
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The parameter a measures a contribution of the second derivative to 
the harmonic oscillator. Therefore, we can expect the classical trajectory to 
behave just like that of the harmonic oscillator when a is small. 

The Euler-Lagrange EOM of the Lagrangian (|4.8p are given by eq. (|2.2p , 

o = Ti-py dL 

= -m 2 q-q-a 2 "q (4.9) 



m 2 + D 2 + a 2 D A )q = (4.10) 



or, equivalently, 



It is not difficult to find clasical solutions to the EOM in eq. (|4.10|) . When 
searching for the classical trajectory in the oscillatory form q c \ = exp(iAi), 
the EOM reads 

^ m 2_ A 2 + a 2 A 4^ el A t =Q (4n) 

and, therefore, we have 



1 ± Vl - 4a 2 



2 



A' (4.12) 

When A is real, the Lagrangian L^PU) is an extension of the harmonic 
oscillator indeed. Hence, we need the condition 

< am< - (4.13) 

It means that the Lagrangina Lpjj has the oscillating solution which is 
similar to the trajectory of the harmonic oscillator. A general solution reads 

q(t) = A + cos(A+t) + B + sin(A+t) + A_ cos(A_t) + 5_ sin(A_t) (4.14) 

where A + , B + , A_ , 5_ are the integration constants, and 



The values of the constants (A + , B + , B_) are determined by the bound- 
ary conditions. 

The Hamiltonian formulation for the Lagrangian (|4.8p can be obtained 
by the Ostrogradski method. The generalized coodinates and momenta are 
given in Sec. 2, ie. 

dL dL 

Qi = q and Pi = — - D— 

dq dq 

dL 

Q 2 = q and P 2 = — (4.16) 
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which imply 



Pi = q + a 2 'q 



P 2 = -a 2 q (4.17) 
The Hamiltonian is given by eq. (|2.12p . ie. 

H = P 1 {DQ l ) + P 2 A-L 

= PiQ2-^Pi-\Ql + ^Ql (4.18) 

or, equivalently, 

H = a 2 qq - y<f + X -q 2 + ^-q 2 (4.19) 

Since the Hamiltonian does not evolve with time, we can find the energy 
by substituting q(t) of eq. (|4.14p at t = into eq. ()4.19f) . as well as q,q,q 
and Q at t = 0, ie. 

9(0) = A + + A^ 
9(0) = B+\+ + B_\_ 

q(0) = -A + \\- A_\ 2 _ (4.20) 
g(0) = -B + X 3 + - £_A 3 

It is now straightforward to calculate the Hamiltonian ()4.19p . We find 

H = a 2 9(0)V(0)-^(0) 2 + ^(0) 2 + ^9(0) 2 (4.21) 
= ^W 1 - ^ 2 m 2 {A\ + B\ ) - ^A 2 \/l - 4a 2 m 2 (A 2 + B 2 _ ) 

To get the Hamiltonian formulation in the harmonic coordinates, we 
begin with the EOM in the form (|4.10p . whose differentioal operator F(D) 
is defined by 

D 2 a 2 D A 

F(D) = 1 + — + ^- (4.22) 

It can be factorized as 

( D 2 



F(D)=[l + ^\n+^) (4.23) 



where X± are given by eq. (|4.15j) . Therefore, the harmonic coodinates are 
given by 



<?+ -(1+^)9 and Q_ = (l + ^ ) q ( 1.2 1 ) 
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The constants rji of eq. (|3.13p can be computed as follows. We have 

dF 



dD 2 



1_ 2a 2 D 2 



so that 



dF 



dD 2 



4 

m 2 



D 2 =-A| 



(1 - 2a 2 \ 2 ±) 



-i 



\2 
. A =fc 

m 2 



\/l — 4a 2 m 2 ^ 



-i 



± 



Al\/1 - 4a 2 m 2 



Therefore, the generalized momenta in eq. (|3.24|) are 



P ± = ^DQ ± 
Xj 



± 



m 



Ai VI - 4a 2 



± 



and the Hamiltonian is given by 

A? 



if 



3 



2m 2 riA 3 



Pf + 



m~ 



E 

3=± 3 



Vj " 2 

4 

m 



{DQjf + XjQj 



(4.25) 



(4.26) 



(4.27) 



(4.28) 



where we have substituted the classical solution (|4.14p . 
The harmonic coodinates (|4.24|) read 

Q+ = A+ fl-^cos(A+f) + B+(l-^±) sin (X+t) (4.29) 

Q_ = A_(l - ^r-) cos(A_i) + 5_^1 - sin(A_f) (4.30) 



where 



A 



4 



±— iyl - 4a 2 m 2 



(4.31) 
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Hence, we find 



Q± = ±-±v / l-Aa 2 m 2 Ia± cos(A±t) + B± sin(A±t) J (4.32) 
Substituting them into the Hamiltonian ()4.28p . we get 

H = ^X 2 + Vl ~ ±a 2 m 2 (A 2 + + B 2 + ) - ^A 2 a/1 - 4a 2 m 2 (A 2 _ + B 2 ) (4.33) 

Equations (|4.22p and (|4.33p are the same. Therefore, we conclude that the 
Hamiltonian formulation by the Ostrogradski method is consistent with the 
Hamiltonian formulation in the harmonic coordinates, as they should. 

The integration constants (A + ,B + ) correspond to the harmonic oscil- 
lator with positive energy, while the integration constants (A- , B- ) corre- 
spond to the harmonic oscillator with negative energy. 



5 Boundary conditions and spectrum 



Going back to the Lagrangian (|4.8p . let us consider its action over a finite 
time period T, 

S[q}= f dtL PU (5.1) 
Jo 

with the trajectory q being a sum of the classical trajectory q c i and the 
displacement q, q = q c i + q. In quantum theory, the displacement q is a 
quantum coordinate. The action can be rewritten as 



q d q-a 2 q c iq+a 2 q d q 



S[q) = S[q d ]+S[q}- J dt (^q d +m 2 q c i+a 2 "q c ^jq+ 

(5-2) 

Here the first term is the action of the classical trajectory q d , and the second 
term is the action of the quantum part q. The integrand of the third term 
vanishes because the classical trajectory is a solution of the (Euler-Lagrange) 
EOM. The fourth term depends on the boundary. However, if the boundary 
condition on q is given by 

J: g(0) = 0, q(T) = 0, g(0) = 0, q(T) = (5.3) 



the fourth term in eq. (|5.2p also vanishes. That boundary condition is the 
same as that of 

^ : ?(0) =g , q{T) = q T , q(0) = %, q(T) = q T (5.4) 
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which was proposed in ref. |21| . The quantum action now takes the form 

%1 . 



■ijf dtq(^D 2 + m 2 + a 2 D 4 ^q+^ 



-] T 



qq — a 2 qq + a 2 qq 



o 

(5.5) 



where the (last) boundary term vanishes due to the boundary condition 

(S3}. 

The boundary term in eq. ()5.5j) also vanishes by another boundary con- 
dition, 

J': g(0) = 0, g(r) = 0, |(0) = 0, if(T) = (5.6) 

As a result, the action (|5.5p takes the Gaussian form, which is quite 
appropriate for a path integral quantization with the Gaussian functional 

dt q (d 2 + m 2 + a 2 DA q (5.7) 



-i 

2 Jo 



Let us now compute the spectrum of the operator D 2 + m 2 + a 2 D 4 . For this 
purpose, we need to find the solutions u k to the eigenvalue equation 

(D 2 + m 2 + a 2 D A )u k (t) = ku k (t) (5.8) 

with the eigenvalues k. A general solution is 

Uk(t) = A\ cos(w + t) + A2 sin(cj + t) + A3 cos(o;_t) + A4 sin(u;_t) 



W± = V 2o^ (5 - 9) 

where A±, A2, A3, A4 is the constants of integration. The function q can be 
expanded in terms of u k , 

dk u k (t) (5.10) 



The spectrum of k is now determined by appying the physical boundary con- 
ditions (|5.3[) or (|5.6p to u k in the form of eq. (|5.9|) . Applying the boundary 
condition ()5.3p at t = yields 

g(0) = A x + A 3 = 0, £(0) = A 2 u + + A 4 w_ = (5.11) 

The boundary condition (|5.3p at i = T then takes the form 

q(T) = Ay cos(cj+T) + A2 sin^+T) — Ai cos(w_T) 

-A 2 — sin(w_T) = 

q(T) = — A\uj + sin(cj+T) + A200+ cos(w + T) 

+Aiw_ sin(w_T) — j42W + cos(w_T) = 
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In particular, the determinant of the matrix on the left side of this equation, 



det 



COsluJ+T) — COS I U!-T) 



\ 



-UJiUJ- 



■uj + sin(w + T) + w_ sin(cj+T) u + 

1 2 



oj- sin(w + r) — uj + sin(w_T) 
cos(w + T) — cos(w_T) 



J/ 



sin 2 (uj i T) +sin 2 (w_T) 



cos(cj + T) — cos(w_T) 

=2o; + o;_ 1 — cos(o; + T) cos(w_T) 

must vanish. We find 

1 — cos(w_|_T) cos(w_T) 

1 — cos(a;_|_r) cos(w_T) 
1 — cos(u;+T) cos(w_T 



— (w+ + u?L ) sin(w + T) sin(a;_T) 
{uj\ + u?L ) sin(w + T) sin(w_T) 



2u> + u>- 
2\/m 2 — 



(5.12) 



+ uPl ) sin(a;_|_r) sin(a;_T) 



a 



■ sinfw+r) sin(a;_T) 



1 



a\/m 2 — k 



sm((jj+T) sin(a;_T) 



(5.13) 



where io±{k) ar given by eq. (|5.9p . Apparently, there is no simple solution 
here. 

When employing the boundary conditions (|5.6p with eq. (|5,9p on it^ , the 
boundary condition in t = yields 



q(0) = A 1 + A 3 = 0, |(0) = -A lW 2 - A 3 oo 2 _ = 



(5.14) 



so that we find A\ = A 3 = when w + 7^ cj_. Now the boundary condition 
at t = T reads 

q(T) = A 2 sin(w + T) + A 4 sin(u;_T) = 

£(T) = -A 2 w 2 sin(u; + r) - A 4 lo 2 _ sin(w_T) = (5.15) 

To get a nontrivial solution, the correspending determinant must vanish, 
which yields the condition 



(w+ — w 2 ) sin(w + T) sin(ci;_r) = 
Since oj + 7^ w_ , we find 

(w + T) =0 or sin(u;_T) = 



sin 



(5.16) 



(5.17) 
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It means 

u + = — or w_ = — (where n is an integer] (5.18) 

and u± are the solutions to the equation 

x 2 + m 2 + a 2 x A = k (5.19) 

Therefore, the spectrum of k with the boundary condition srf 1 has the simple 
form 

k= (^y +m * +a ^y (5.20) 

6 Canonical quantization and instabilities 

In this section we recall about istabilities and ghosts in the quantum PU 
oscillator [13] . The most straightward way is based on identifying the energy 
rasing and lowering operators [H]. The classical solution (|4.14p can be 
rewritten to the form 

q(t) = ±(A + -iB + y x + t + ±(A + + iB + )e- a + t 

- iBJ)e iX - 1 + \{A-+ iBJje^* (6.1) 

Since the A_ modes have negative energy, the lowering operator must be 
proportional to the (^4_ — iBJ) amplitude. Similarly, since the A+ modes 
have negative energy, the raising operator must be proportional to the {A + + 
iB + ) amplitude, ie. 

a± ~ A± db iB± 

~ ^(1 ±y/l- 4a 2 m 2 )Qi ± iP\ T |(1 4a 2 m 2 ) - \±P 2 

(6.2) 



where we have used 



and 



A ± _ (6 . 3 ) 



as 



well as0 



5± = A ± (A 2 -A 2 ± ) (6 - 4) 



Qi = g (6.5) 



lr The canonical variables were calculated at the initial time value because the operators 
in Schrodinger picture do not depend upon time. 
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Qi = go (6.6) 

P 1= q + a 2 q (6.7) 

P2 = -a 2 q (6.8) 
It is now straightfoward to derive the commutation relations, 

[a±,c4] = l (6.9) 



The next step depends upon physical interpretation [14j . 

(I) The 'empty' (or 'ground') state may be defined by the condition 



a+ In) = al \n 







(6.10) 



Then the 'empty' state wave function £l{Qi,Q2) (in the Q-representation, 
with P = —id/dQ) reads 



0(Q 1) Q 2 ) = Nexp 



imaQiQ 2 



(6.11) 



and is infinite or not normalizable, because the size of the wave function 
gets bigger with the increase of Q2 , so that the integral over the whole space 
diverges. 

In addition, when the eigenstate jV+,iV_) with the eigenvalues N = 
(N + ,N-) is defined by 



,t 



\N+,N- 



the norm of the (0, 1) state is given by 

<0,1|0,1> = (n\ala-\n) 



(6.12) 



(-l + a_at) IQ; 
- <tt\&> 
-1 



(6.13) 



which is a ghost. The non-normalizable quantum 'states' are physically 
unacceptable, so the interpretation (I) should be dismissed |14j . 

(II) It is, however, possible to treat all particles (with positive or negative 
energy) as the truly ones by defining the 'empty' state differently, namely, 
as 

a±|fi) = (6.14) 

In this interpretation the negative energy can arbitrarily decrease and the 
Hamiltian is unbounded from below. The 'empty' state solution £l(Qi,Q2) 
in the Q representation is now given by 



n(Qi,Q 2 ) = Nexp 



vT — 4a 2 m 2 
' 2(A_ + A+) 



(X+X-Qj + Ql) + imaQxQi 



(6.15) 
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and is apparently finite or normalizable, because the first term in the expo- 
nential is negative. 

The eigenstate \N+, N-) of the eigenvalues N = (N+, iV_) is now given 

by 

, 

\N+,N-) = —^=—^=\n) 6.16 

^/n^.^/nJ. 

while the norm of the (0, 1) state is 

<0,1|0,1> = (n\a-al\n) 

= {n\ (1 - at a_) \n) 

= <n\n> 

= 1 (6.17) 

ie. it is not a ghost. 

In the correct physical interpretation (II) the correspondence principle 
between the classical and quantum states is preserved, but the system has 
indefinite energy. When interactions are switched on, mixing the negative 
and positive energy states would lead to instabilities in the classical theory, 
and the exponentially growing and decaying states in quantum theory [26, 
127] . Excluding the negative energy states would lead to the loss of unitarity 

m\. 



7 Path integral quantization and Forman theorem 

The idea of ref. |21j is to define the quantum theory of the PU oscillator 
as the Euclidean path integral and then Wick rotate it back to Minkowski 
case. It makes sense since the Euclidean action of the PU oscillator — see 
eq. (|8.3p below — is positively definite. It can also make the difference to 
the canonical quantization and the Ostrogradski method (Sec. 2) when one 
integrates over the path only, but not over its derivatives. 

Let us first recall some basic facts about a path integral in QFT, accord- 
ing to the standard textbooks in Quantum Field Theory - see, for example, 
ref. [2H]. 

The definition of the probability amplitude for a one-dimensional quan- 
tum particle by Feynman path integral is given by 



rib 






Z(qb,t b ;q a ,t a ) = & 
J q a 


gexp 


~ ^dtL 



where the integration goes over all paths q(t) between q a and After Wick 
rotation 

t->t = -ir (7.2) 
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the path integral takes the fori 



Z{qb,tb;q a ,t a ) 



it 



cItLe 



(7.3) 



It is called the Euclidean path integral. In the case of the PU oscillator the 
Euclidean path integral is Gaussian. Let us recall some basic properties of 
the Gaussian integrals. 

The simplest Gaussian integral reads 



dxe 



a > 



It can be easily extended to a quadratic form in the exponential as 



dxe 



-ax —bx 



- exp fy 



(7.4) 



(7.5) 



It can also be easily extended to the case of several variables with the 
diagonal quadratic form as 



J [d n x] exp ( — aixj J = r 



(7.6) 



where we have introduced the normalized measure [dx] = dxj 'y/w. 

By diagonalizing a generic (non-degenerate) quadratic form, one can 
prove a general finite-dimensional formula, 

/ [d n x] exp (-s* Ax - tfx) = 1 exp ( ^A'H ) (7.7) 

Joo lli=l *i V 4 / 

= exp (\\fA- x b 

detA 1 \A 

Finally, when formally sending the number of integrations to infinity, 
one gets the Gaussian path integral, 



rib 




rtb 


/ ^ 


Iqexp 








. Jt a 



dt{q{t)F{D)q{t)+q{t)J{t)) 



exp 



yjDetF(D) 

where DetF(D) is now the functional determinant 



dtJ(t)F~ 1 (D)J{t) 



(7.8) 



The sign factor in the Wick rotation is chosen to make the path integral converging. 
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A generic functional determinant diverges since it is denned as the prod- 
uct of all the eigenvalues in the spectrum of a differential operator. There- 
fore, one needs a regularization. It is most convenient to use the zeta func- 
tion regularization in our case — see, for example, ref. [29J for a comprehen- 
sive account. The Riemann zeta function is defined by 



oo ^ 



(7.9) 



71=1 



in the convergence area of the series. It is then expanded for Re(s) > 1 by 
analytic continuation. It is often useful to employ an integral representation 
of the zeta function in the form 



-I POD x 



Tit 



(7.10) 



where the (Euler) gamma function has been introduced, 

r(s) = 



dtt s - l e- 1 



(7.11) 



Equation (|7.1U|) allows one to define the zeta function for an elliptic 
operator L as 



as\L) 



where tre tL is given by 



tre 



-tL 



tr 



r(s) Jo 



dtt s -hre~ tL 



(7.12) 



-\ 2 t 



71=1 



-X n t 



(7.13) 



in terms of the positive eigenvalues X n of L. One easily finds 

oo 1 



n=l A « 



n=0 



-s hi A„ 



(7.14) 



Differentiating both sides of this equation with respect to s at s = 0, one 
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finds 



dC(s\L) 



oo 



ds 




■n 



s=0 



n=l 



oo 



In H A 



■n 



n=l 



In DetL 



(7.15) 



so that the functional determinant of an elliptic operator L is given by 



The zeta function regularization of the right hand side of this equation is 



where we have introduced the regularization parameter e and the dimension 
parameter /x. 

The zeta-function renormalization amounts to deleting the first term in 
eq. (|7.17p . since it UV-diverges in the limit e — > 0, as well as the third term 
since it IR-diverges in the limit [i — > 0. 

To put equation (|7.17p into a more explicit form, without resorting to 
the spectrum of the differential operator, it is convenient to use Forman's 
theorem [23] :| 

Let and are the differential operators defined by 



DetL = e -C'U=o 



(7.16) 




e 

i(l + 6ln/, 2 )(C(0|L) + eC / (0|L)) + O(6 2 ) 
-C(0|L) - C'(0|L) - ln/i 2 C(0|L) + 0(e 2 ) (7.17) 




(7.18) 



3 Forman theorem is an extension of the GePfand-Yaglom theorem. 
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over the domain [0, T\. Consider a linear differential equation 

Kh{r) = 

with a boundary condition 



(7.19) 



: M 



( Ho) \ 

hW{0) 

\ h^-^io) J 



+ N 



( HQ) \ 

hW{0) 

\ h^-^io) J 



(7.20) 



and take the boundary condition to be smoothly connected to g/. The 
time evolution operator Yk{t) is introduced as 



(7.21) 



/ Kt) \ 


= Y K (r) 


/ mo) \ 






V ^"-^(o) / 



so that the boundary condition can be written to 

/ HO) 
(M + NY K (T)) : 

\ h( n ~ l X0) 

The Forman theorem is given by the statement: 

DetK^ _ det(M + NY K (T)) 
DetKj ~ det [M + NY R (T)) 



(7.22) 



(7.23) 



This theorem is effective for finding the functional determinant of the 
operator K with unknown spectrum by connecting it to the one with a 
simple spectrum via changing the boundary conditions. 

8 Path Integral of PU Oscillator 

The Euclidean path integral of the PU oscillator over a domain [0, T] was 
calculated in refs. [21 1 122 ^ 125]. Here we confirm the results of ref. |22| by our 
calculation. 

The path integral of PU oscillator with the action 



2 ^{t?-\m 2 



after the Wick rotation (t —> it) takes the form 

Z(q T ,T;q ,0)= / @qexp(-S E ) 
J an 



1.2) 
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where the Euclidean PU action is given by 

; / 1 . , , 9 m 2 , N 9 a 2 



S E = J dt ^-q(tf + ~YQ(t) 2 + Y q(t) 2 ) (8-3) 

This Se is positively definite, so that the Euclidean path integral is well 
defined. 

Since our discussion of the classical theory (Sec. 4), the integral tra- 
jectory is a sum of a classical trajectory q c \ and quantum fluctuations q, 
Q = Qci + Q- Accordingly, the action can be also written down as a sum, 

S E [q}=S cl + S[q] (8.4) 

and the path integral of the PU oscillator takes the form 

r° 

Z(q T ,T;q ,0) = e- S <* 3>qexp (S[q\) (8.5) 
Jo 

where the quantum action S[q] is given by 

Sm = \[dtq{a^-^ 2+ rn 2 )q (8.6) 

after integration by parts. 

Let us denote the differential operator o?-^ — ^7+ m - 2 with the boundary 
condition as K^. Then the path integral can be written down in the form 

Z(q T ,T;q ,0) = e- S ^ J ®qexp J dtqK^q\ (8.7) 

The path integral of the PU oscillator is Gaussian and, therefore, can be 
computed along the lines of Sec. 7 as 

N 

Z(q T , T; g , 0) = exp (S d ) (8.8) 

where N is the normalization constant. The classical part Sd was found 
in ref. [21], and it is quite involved. The functional determinant is the key 
part of a quantum propagator of PU oscillator, which is of primary physical 
interest. It can be computed by the use of Forman theorem (Sec. 7). 
First, one calculates the time evolution operator Yk- It is given by 



Y K (t) 



( ui(t) u 2 (t) u 3 (t) u 4 (i) \ 

U\{t) u 2 {t) u 3 (t) u 4 (t) 

u'i(i) u 2 (t) u 3 (t) u 4 (t) 

V ui(t) ii 2 (t) u 3 (t) ul{t) J 



(8.9) 



where 

K Ui {t) = (i = 1,... ,4) (8.10) 
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and the inital condition is 



n(o) 



(8.11) 



(8.12) 



The operator is equal to K^, so they have Yx(t) is common. 
By solving the equation Ku-i = for Ui with 

2 d 4 d 2 2 

one gets its general solution in the form 
Ui(t) = Ai sinh(A+t) + A cosh(A+i) + C L sinh(A_t) + A cosh(A_t) (8.13) 
The boundary condition Yrt(O) = 1 amounts to the relations 

A + A = S u 
X + Ai + A_Cj = 5 2 i 



X\Bi + A 2 A = 5 3i 
X\Ai + \^_Ci = 5a 

Therefore, the solutions are 

A 2 _ . , A 2 



(8.14) 
(8.15) 
(8.16) 
(8.17) 



Ul 



U2 



U3 



U4 



A 2 



cosh(A + t) + 



A 2 - A 2 



cosh(A_t) 



A + (A 2 _-A 2 ) 
1 



sinh(A + i) + 



A 2 



A 2 



1 



cosh(A + i) 



A 2 (A 2 h 
1 

A 2 - A 2 . 



A + (A 2 _ 



a 2 + : 



sinh(A + t) + 



"A 2 , 
cosh(A_t) 
1 



sinh(A_t) 



(8.18) 
(8.19) 
(8.20) 



—2 ^smh(A_t) (8.21) 

A 2 (A + 



A 2 J 



Next, one writes down the boundary conditions si and si in terms of 
the matrices M and N appearing in the Forman theorem. The boundary 
condition si is 



si : q(0) = , q(T) = , g(0) = , q(T) = 
so that its matrices M and N are given by 



M 



(8.22) 



and 



N = 



( 1 








o \ 





1 




















V o 








o / 


( 








o \ 














1 











V o 


1 





o / 



(8.23) 



(8.24) 
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In the same way, the boundary condition s$ is 

d : $(0) = , q(T) = , |(0) = , q\T) = 
so that its matrices M and N are given by 



M 



.25) 



and 



N 



( 1 
















-1 

















V o 








o / 


/ o 








°\ 














1 











\o 





1 


o / 



.26) 



.27) 



Having found M, N and Yk, as well as M, N and Y R , we calculate 



det(M + NY K (T)) 



a 
m 



1 



1 + 2ma 



■ sinh 



VT + 2ma, 
2a 



1 



1 — 2ma 



■ sinh 



\/T— 2ma, 
2a 



3.28) 



and 



det(M + NY R {T)) 
sinh 2 



a 
m 



y/l + 2ma r 
2a " 



sinh 



yjl — 2ma r 
2a " 



.29) 



A calculation of DetK^ goes along the standard lines [211 221 [23], 

oo 

DetK^ = Y[k n 



n=l 
oo 



n=l v 7 



3.30) 



2^2 



,iah 2 f V1 + 2ma r 

2a 



sinh 2 ( Vl T 2ma T 



2a 



By using the Forman formula, we get the final answer 
cfei (M + NYk{T)) 



DetK 



det (M + NY R (T)) 



DetK r 



a 



m 



T 2 



[1 + 2amy 1 sinh 2 



y/l + 2ma r 



2a 



3.31) 



-(1 — 2am) 1 sinh 2 



\/l — 2ma 
2a 



T 
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in full agreement with ref. [22] in its last (v2) version. In the large T>1 
limit one finds 



DetK^ = (8.32) 



a 
m 



1 exp (a/1 + 2ma^) 1 exp (\/l — 2ma^) 



l + 2am (2T) 2 l-2am (2T) 2 

and in the small T limit one gets 

T 2 

Deiif^ « — + 0(T 4 ) (8.33) 

The transition amplitude (or the quantum Euclidean propagator) of PU 
oscillator is given by 



/ 2-7T 

< Qt, Qt; t = T\q , q ;r = >= J DetR ^ ex P [SeIQci]) (8.34) 

The classical Euclidean action Se[<1cI\ was calculated in Appendix of ref. [2l~] . 
It is finite for large T> 1 and behaves like for small T -C 1. Hence, the 
transition amplitude is exponentially suppressed both for small and large T, 
ie. it is normalizable, and the Euclidean path integral is well defined indeed. 



9 Conclusion 

The procedure of calculating Euclidean transition probabilities (for observ- 
ables) in the quantum PU theory was outlined in ref. |21| . The probabilities 
in the Minkowski space can be obtained by analytic continuation. It is, 
therefore, possible to make physical sense out of the quantum PU theory. 

In classical PU theory with interactions, even at a very small value of 
the parameter a > 0, one gets runaway production of states with negative 
and positive energy. However, as was suggested in ref. [21] . the Euclidean 
formulation of the quantum theory implicitly imposes certain restrictions 
that can remove classical instabilities. The price of removing the instabili- 
ties is given by an apparent violation of unitarity [21J. Indeed, integrating 
over the basic trajectory, and not over its derivatives in the Euclidean path 
integral formulation of the quantum PU oscillator given above is not in line 
with the canonical quantization and the Ostrogradski method. By doing it, 
one looses some information and, hence, one loses unitarity. As was argued 
in ref. [21], one can, nevertheless, never produce a negative norm state or 
get a negative probability, so that the departure from unitarity may be very 
small at the low energies (say, in the present universe) , but important at the 
very high energies (say, in the early universe). Of course, it is debateable 
whether the 'price' of loosing unitarity is too high or not. 
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Apparently, the f(R) gravity theories are special in the sense that for 
each of them there exist the classically equivalent scalar-tensor field theory 
without higher derivatives, under the physical stability conditions. Still, as 
the quantum field theories, they may be different. Hence, it may be possible 
to quantize f(R) gravity without loosing unitarity. Figuring out the details 
is still a challenge. 
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